Jet tomography of high-energy nucleus-nucleus collisions at next-to-leading order 
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We demonstrate that jet observables are highly sensitive to the characteristics of the vacuum 
and the in-medium QCD parton showers and propose techniques that exploit this sensitivity to 
constrain the mechanism of quark and gluon energy loss in strongly-interacting plasmas. As a first 
example, we calculate the inclusive jet cross section in high-energy nucleus-nucleus collisions to 
0{al). Theoretical predictions for the medium-induced jet broadening and the suppression of the 
jet production rate due to cold and hot nuclear matter effects in Au+Au and Cu+Cu reactions at 
RHIC are presented. 
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PACS numbers: 13.87.-a; 12.38.Mfi; 25.75.-q 

Hadronic jets pj], collimated showers of energetic final- 
state particles, have long been regarded as the main tool 
to understand hard scattering (Q^ ^qcd) processes 



uncertainty in the extraction of the plasma properties, 
such as its density and transport coefficients [8|. 

To take full advantage of jet physics in reactions with 
in e+ -I- e~, semi- inclusive deeply inelastic scattering, and ultra-relativistic nuclei, calculations at 0{al) are re- 



hadron-hadron collisions from first principles in Quantum 
Chromodynamics (QCD). Due to their large production 
rate [H , jets are easily accessible by experiment and their 
landmark discovery has, arguably, stimulated some of the 
most important developments in the perturbation theory 
of strong interactions. At present, cross sections for pro- 
cesses involving jets are routinely calculated at next-to- 
leading order (NLO) and next-to-next-to-leading order 
results are also becoming available Q- This remarkable 
theoretical accuracy has made precision jet studies the 
optimal method for finding new physics beyond of the 
Standard Model at very high energies 0\. 

An unprecedented possibility exists today to extend 
the theory of jets to energetic reactions with large nuclei 
(A-l-A) [J]. Experimental advances at this interface be- 
tween particle and nuclear physics have already allowed 
proof-of-principle measurements of jets at the Relativis- 
tic Heavy Ion Collider (RHIC) Q. In this Letter we 
present the first calculation of jet cross sections and jet 
sub-structure for Au-}-Au and Cu-|-Cu reactions at RHIC 
that includes both NLO perturbative effects and effects 
of the nuclear medium. We demonstrate how jet observ- 
ables can be used to gain insight into the mechanisms of 
parton interaction and energy loss [6| in hot, dense, and 
strongly-interacting quark-gluon plasmas (QGPs) that 
are created in such reactions. We determine quantita- 
tively the relation between the in-medium modification 
and broadening of parton showers in the QGP and the 
suppression in the observed cross section as a function of 
the jet cone radius R. Owing to their differential nature, 
jet observables will soon place stringent constraints [E, 0| 
on the approximations and theoretical model assump- 
tions that underlay the current competing approaches to 
quark and gluon energy loss in hot and dense QCD mat- 
ter. Such discriminating power cannot be achieved with 
measurements of leading particle quenching. This will, in 
turn, help eliminate the uncomfortably large systematic 



quired In fact, the lowest order QGP-induced par- 
ton splitting is also manifested in experimental observ- 
ables at the same 0{al). With the shortcut notation 
d{ET,y,(j)}n = Y{i=i^n'^yiY{j=2,n^Tjd(j)j for the inde- 
pendent final-state parton variables in the collinear fac- 
torization approach, the simplest inclusive jet cross sec- 
tion at NLO can be expressed as follows Q: 
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Here, ETi,yi,(f>i are the transverse energy, rapidity, and 
azimuthal angle of the i-th particle {i = 1,2,3), respec- 
tively, and (j[2 — >■ 2], a[2 — 3] represent the produc- 
tion cross sections with two and three final-state par- 
tons. In Eq. ([U S'2, S'a are phase space constraints and 

^2 = ELi = ELi HEt, - ET)S{y^ - y) identifies 
the jet with its parent parton. Hence, only at NLO can 
the dependence of the experimental observables on the 
jet cone radius R and the jet finding algorithms 0, 
be theoretically investigated. For an angular separation 
Rij — \J (j/i — yjY' + ~ 4>iY' ^ defined for any possible 
parton pair (i, j), 
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In Eq. (HI) i?,c = min (^RsepR, ^,^^Et,,e^ 
when two partons should be recombined in a jet. Here, 
1 < Rsep < 2 is introduced to take into account features 
of experimental cone algorithms, employed to improve 
infrared safety. Eq. ([2]) establishes a correspondence be- 
tween the commonly used jet finders and the perturbative 
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FIG. 1: Comparison of the NLO calculation to STAR exper- 
imental data 11] on the inclusive jet cross section for R=0.4 
(top panel). The variation of the jet cross section with the 
cone size R for Et ~ 10, 40 GeV around midrapidity at RHIC 
is also shown (bottom panel). 



calculations to ©(af) with the goal of providing accu- 
rate predictions for comparison to data. For example, 
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2 yields a midpoint cone algorithm and R 



Sep 



corresponds to the kx algorithm [2|, [10| . R is the cone 
size or parton separation parameter, respectively. 

In the top panel of Fig. [1] we show a comparison of 
the NLO calculation Q of the inclusive jet cross sec- 
tion at y/s — 200 GeV p+p collisions at RHIC to 
the STAR experimental measurement which uses a mid- 
point cone algorithm [11] in the pseudorapidity range 
0.2 < rj < 0.8. Very good agreement between data and 
theory is achieved with a standard choice for the renor- 
malization and factorization scales /i^ = A*/ — ^t- Vari- 
ation of these scales within {Et/'2, 2Et) leads to less than 
(-1-10%, —20%) variation of the jet cross section. The 
bottom panel of Fig. [T] illustrates the significant depen- 
dence of da^'^'^ /dydEr on the cone size R, which, even in 
p-|-p reactions, can exceed a factor of two. Analytically, 
the ln(i?/i?o) scaling of the cross section can be under- 
stood from the 1 /r angular behavior of the perturbative 
QCD splitting kernel. Further insight into the underly- 
ing parton dynamics can be gained by examining the jet 
sub-structure, often characterized by the differential jet 
shape: 
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evaluated as follows 0, '■ 

r^'-{r, R) = ^con(r-, R) (Psudakov(r, i?) - 1) + ^Lo(r, R) 
+V'»,LO {r, R) + ?Apc (r, R) + V'»,pc {r, R) . (4) 

In Eq. (j4l) the first term represents the contribution from 
the Sudakov-resummed small-angle parton splitting; the 
second and third terms give the leading-order final- 
state and initial-state contributions, respectively; the 
last two terms come from power corrections cx Qo/Et, 
Qo — 2 — 3 GeV, when one integrates over the Landau 
pole in the modified leading logarithmic approximation 
(MLLA). This approach was shown to provide a very 
good description of the differential intra-jet energy flow 
at the Tevatron as measured by CDF II [l2]. Thus, 
reliable predictions for the jet sub-structure in p-l-p re- 
actions at RHIC can be obtained and used as a base- 
line to study the distortion of jet shapes in more com- 
plex systems, such as p-|-A and A-l-A. An example of 
■0™'= (r,i? ^ 0.4) for a quark jet of Et = 30 GeV is 
shown in Fig. [2] 

When compared to a parton shower in the vacuum, the 
medium-induced quark and gluon splittings have notice- 
ably different angular and lightcone momentum fraction 
dependencies [i, 13|. In particular, for energetic partons 
propagating in hot and dense QCD matter, the origin of 
the coherent suppression of their radiative energy loss, 
known as the Landau-Pomeranchuk-Migdal effect, can 
be traced to the cancellation of the collinear radiation at 
r < mD/ {ijj{rnDT^g,ET)) flS*]. Here, the Debye screen- 
ing scale m_D = gT^\ -I- Af//6 and (i^) — few GeV. 
Thus, the medium-induced component of the jet, which 
is given by the properly normalized gluon bremsstrahlung 
intensity spectrum '\\}™'^^(r^T() oc dF^^ j dujdr within the 
cone, has a characteristic large-angle distribution away 
from the jet axis. This is illustrated in Fig. [2] for cen- 
tral Au-|-Au and central Cu-KCu collisions at RHIC. We 
emphasize that accurate numerical simulations, taking 
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FIG. 2: The differential jet shape in vacuum 'i/)™'^ (r, _R) is 
contrasted to the medium-induced contribution t/i™'^'*' (r, _R) 
by a Et = 30 GeV quark in Au-|-Au and Cu-|-Cu collisions 
at y^sjvjv = 200 GeV. The insert illustrates a method for 



In Eq. ([2]) i stands for the sum over all particles in this studying the characteristics of these parton showers, 
jet and /g '0(r, R}jdr — 1. Analytically, jet shapes can be 
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into account the geometry of the heavy ion reaction, the 
longitudinal Bjorken expansion of the QGP, and the con- 
straints imposed by its experimentally measured entropy 
density per unit rapidity [J] , have been performed for all 
physics results quoted in this Letter. 

One can exploit the differences between the vacuum 
and the in-medium parton showers by varying the cone 
radius R (i?i,jet < R) and a cut p??'" [Et, > for the 

particles "z" that constitute the jet, to gain sensitivity 
to the properties of the QGP and of the mechanisms of 
parton energy loss in hot and dense QCD matter. This is 
illustrated in the insert of Fig. [2j The most easily accessi- 
ble experimental feature of jet production in nuclear col- 
lisions is, arguably, the suppression of the inclusive cross 
section in heavy ion reaction compared to the binary col- 
lision scaled, oc (-/Vbin), production rate in elementary 
nucleon-nucleon reactions 141: 



V^^bin; dyd-^ET 



(5) 



Eq. ([U defines a two dimensional jet attenuation pattern 
versus R and p™'" for every fixed Et- In contrast, for the 
same Et , inclusive particle quenching is represented by a 
single value related to the R ^ and p™'" ^ (lu) limit in 
Eq. ([5]). Thus, jet observables are much more differential 
and, hence, immensely more powerful than leading par- 
ticles and leading particle correlations in their ability to 
discriminate between the competing physics mechanisms 
of quark and gluon energy loss in dense QCD matter and 
between theoretical model approximations to parton dy- 
namics in the QGP. 

We now focus on the first complete theoretical result 
at NLO for R'2li versus the jet cone size R for Au+Au 
and Cu+Cu collisions at RHIC. We use the reaction op- 
erator approach to non-Abelian energy loss [l4|, in the 
limit of weak coupling between the jet and the plasma 
with as ^ 0.3, to evaluate the probability distribution 
Pq^g{t,E) that quarks and gluons, respectively, will lose 
a fraction of their energy e = iOijE due to medium- 
induced bremsstrahlung. Next, we determine the fraction 
of this energy that will be redistributed inside the jet: 



duj 



dPJ 



dujdr 



Jo Jo duidr 



(6) 



While such redistribution may affect the jet shape, it 
will not affect the jet cross section. For example, when 
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1) final-state QGP- 



induced effects to inclusive or tagged jet cross sections 
vanish. Parton interactions in the strongly-interacting 
plasma, however, are not the only many-body QCD ef- 
fects that will alter the measured jet cross section. Cold 
nuclear matter (CNM) effects prior to the QGP forma- 
tion [15] must also be included in accurate theoretical 
calculations of hard probes production in nuclear colli- 

We find that in 



sions and we first evaluate 
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FIG. 3: Transverse energy dependent nuclear modification 
factor R'2\ for different cone radii i? in 6 = 3 fm Au+Au 
(top panel) and Cu+Cu (bottom panel) collisions at ^/sNN = 
200 GeV. Inserts show ratios of jet cross sections for different 
R in nuclear reactions versus Et- 



the kinematic region of interest, 10 GeV < Et < 50 GeV 
around midrapidity at RHIC y^SNN — 200 GeV colli- 
sions, the EMC effect and initial-state energy loss [l3| 
play a dominant role. Next, we determine the relative 
fractions Uq^g of quark and gluon jets in this inclusive 
cross section (?Xo + rig — 1). These are well defined at 



da]- 



leading order [16| and separation of the inclusive cross 
is necessary to properly describe 



section into 



d'-'Erdy 



parton energy loss in the QGP, which scales with the 
quadratic Casimir in the corresponding representation of 
SU(3) (Ca/Cf = 2.25). At NLO there exists an ambi- 
guity of 0{as) in this separation which has a very 
small effect on inclusive jet observables. 

We calculate the medium-modified jet cross section per 
binary nucleon-nucleon scattering as follows (p™'" = 0) : 



1 da^^jR) 
(iVbin) d^ETdy 



f de J2PqJe,E) 



j^CNM,NLO 
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In Eq. ([7]) (1 — fq^g) ■ e represents the fraction of the en- 
ergy of the parent parton that the medium re-distributes 
outside of the cone of radius R. The measured cross 
section is then a probabilistic superposition of the cross 
sections of protojets of initially larger energy Elj. — 
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{r/R) 


Vacuum 


Medium 


Total 


A 


Au+Au 


0.271 


0.601 


0.283 


4% 


Cu+Cu 


0.271 


0.640 


0.272 


0.4% 



TABLE I: Mean relative jet radii (r/R) in the vacuum, 
for complete parton energy loss in the medium, and for the 
realistic case of jets in Au+Au and Cu+Cu collisions at 
RHIC. We considered a radius R = 0.4 and transverse en- 
ergy Et = 30 GeV at y/sJ^N = 200 GeV. The fractional 
QGP-induced broadening A{r/R) is also shown. 



Et/{1 — (1 — fq,g) ■ e). Our results for the nuclear modifi- 
cation factor of inclusive jets R'2\ central Au+Au and 
Cu+Cu collisions with ^/snn = 200 GeV at RHIC are 
presented in Fig. [3] Each band illustrates a calculation 
for a 20% increase in the rate of parton energy loss 
(lower bound) relative to our default simulation (upper 
bound) . Experimental data on leading tt*' suppression for 
these reactions is only included for reference. A contin- 
uous variation of R^2\ with the cone radius R is clearly 
observed in Fig. [3] and shows the sensitivity of the inclu- 
sive jet cross section in high-energy nuclear collisions to 
the characteristics of QGP-induced parton shower. For 
R < 0.2 the quenching of jets approximates the already 
observed suppression in the production rate of inclusive 
high-pT particles. It should be noted that in our theo- 
retical calculation CNM effects contribute close to 1 /2 of 
the observed attenuation for Et > 30 GeV. These can be 
drastically reduced at all Et by taking the ratio of two 
differential cross section measurements for different cone 
radii i?i and i?2. A few selected examples are shown in 
the inserts of Fig. [3] 

Inclusive jet cross sections and jet shapes in nuclear 
collisions are closely related Q: 
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should be noted that vacuum and medium-induced 
parton showers become more coUimated with increas- 
ing E^ and the mean relative jet width {r/R) — 
d{r/R){r/R)iP{r/R) is reduced Q. Consequently, 
the striking suppression pattern for jets, shown in 
Fig. ([HI, can be accompanied by a very modest growth in 
the observed {r/R) tot.- We show in Table U the relative 
widths in the vacuum, for a hypothetical case of com- 
plete parton energy loss {Pq,g{i) = 5{\ — e)) that falls 
inside of R {fq^g = 1), and for our realistic simulation 
of Et = 30 GeV jets of i? = 0.4 in central Au+Au and 
Cu+Cu collisions at RHIC. We find that, on average, jet 
broadening A(r/i?) = ((r/i?)tot. - (r/i?)vac.)/(r/-R)vac. 
is < 5%. Larger effects are expected near the core r — >• 



(8) 



It 



and the periphery r — ^ i? of the jet [J|. 

In summary, the principle goal of this Letter is to 
bring into focus the immense possibilities that jets of- 
fer as tomographic probes of the QGP created in ultra- 
relativistic nuclear collisions. To demonstrate the un- 
precedented sensitivity of jet observables to the charac- 
teristics of the vacuum and the medium-induced parton 
showers, we presented first results for the related cross 
sections and shapes as a function of the cone radius R 
at next-to-leading order 0{al) in p+p, central Au+Au, 
and central Cu+Cu collisions at RHIC. Our theoretical 
predictions include a most detailed account of cold and 
hot nuclear matter effects on jet production and distor- 
tion in heavy ion reactions and are of immediate relevance 
to upcoming PHENIX and STAR experimental measure- 
ments. We fully expect that this work will inspire future 
studies of inclusive and tagged jets, jet sub-structure and 
event shape observables. In their entirety, such studies 
will provide first-principles insights into the many-body 
QCD parton dynamics at ultra-relativistic energies and 
shed light on the relative importance of collisional ver- 
sus radiative energy loss and on the applicability of weak 
versus strong jet-medium coupling regimes. 

Acknovifledgments: We thank D. E. Soper and 
F. I. Olness for helpful discussions. 



[2: 



[5: 

[6: 

[7: 

[s: 
[9: 

[10 

[11 

[12 

[13 
[14 

[15 
[16 



G. Sterman, S. Weinberg, Phys. Rev. Lett. 39, 1436 
(1977); R. P. Feynman, R. D. Field, G. C. Fox, Phys. 
Rev. D 18, 3320 (1978). 

J. M. CampbeU, J. W. Huston and W. J. Stirling, Rept. 
Prog. Phys. 70, 89 (2007); S. D. EUis et al. Prog. Part. 
Nucl. Phys. 60, 484 (2008). 

F. I. Olness and D. E. Soper, arXiv:0907.5052 [hep-ph]. 
I. Vitev, S. Wicks and B. W. Zhang, JHEP 0811, 093 

(2008) ; 1. Vitev, B. W. Zhang and S. Wicks, Eur. Phys. 
J. C 62, 139 (2009). 

S. Salur, Nucl. Phys. A 830, 139C (2009); M. Ploskon, 
Nucl. Phys. A 830, 255C (2009). 

X. N. Wang, M. Gyulassy, Phys. Rev. Lett. 68, 1480 
(1992). 

N. Gran, Nucl. Phys. A 830, 797C (2009); D. d'Enterria 

et al.,3. Phys. G 34 (2007) 2307. 

S. A. Bass et al. Phys. Rev. C 79, 024901 (2009). 

S. D. Ellis, Z. Kunszt, D. E. Soper, Phys. Rev. Lett. 64, 

2121 (1990); Z. Kunszt and D. E. Soper, Phys. Rev. D 

46, 192 (1992). 

M. H. Seymour, Nucl. Phys. B 513, 269 (1998). 

B. I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 

97, 252001 (2006). 

D. E. Acosta et al. [CDF Collaboration], Phys. Rev. D 
71, 112002 (2005). 

I. Vitev, Phys. Lett. B 630, 78 (2005). 

M. Gyulassy, P. Levai, I. Vitev, Phys. Rev. Lett. 85, 5535 

(2000); I. Vitev, Phys. Rev. C 75, 064906 (2007). 

I. Vitev, B. W. Zhang, Phys. Lett. B 669, 337 (2008); 

R. Sharma, L Vitev, B. Zhang, Phys. Rev. C 80, 054902 

(2009) . 

M. Dasgupta, L. Magnea, G. Salam, JHEP 0802, 055 
(2008). 



